Order parameter for two-dimensional critical systems with boundaries 

Ivica Res and Joseph P. Straley 
Department of Physics and Astronomy, University of Kentucky, Lexington, KY 40506-0055 



o 
o 
o 

(N 

X) 
CD 

in 



o 

e 

I 

-(— > 



o 

> 

O 

On 
OS 
-I— > 



I 

O 

o 



X 



Conformal transformations can be used to obtain the or- 
der parameter for two-dimensional systems at criticality in 
finite geometries with fixed boundary conditions on a con- 
nected boundary. To the known examples of this class (such 
as the disk and the infinite strip) we contribute the case of a 
rectangle. We show that the order parameter profile for sim- 
ply connected boundaries can be represented as a universal 
function (independent of the criticality model) raised to the 
power ir;. The universal function can be determined from the 
Gaussian model or equivalently a problem in two-dimensional 
electrostatics. We show that fitting the order parameter pro- 
file to the theoretical form gives an accurate route to the de- 
termination of r]. We perform numerical simulations for the 
Ising model and percolation for comparison with these ana- 
lytic predictions, and apply this approach to the study of the 
planar rotor model. 

PACS numbers: 64.60.Cn, 05.70. Jfc, 64.60.Fr 



I. INTRODUCTION 

A system undergoing a continuous phase transition is 
scale invariant at its critical point. An interesting and 
productive generalization of this idea for systems also 
having translational and rotational invariance is to as- 
sume conformal invariance [y^]. Among other things, 
this implies that the magnetization (the one-point cor- 
relation function) for a two-dimensional critical system 
transforms under a conformal transformation w = w{z) 
in the following way: 



< cp{w,w) >= {^)-'-A-'^ < 0(z,z) > 



dz 



dz 



(1) 



Here <> means the statistical average, z — x + iy is 
a complex-valued coordinate, and it is assumed that 
< (f)(z, z) > is nonzero due to a breaking of the (f) field 
symmetry at a boundary. For a spinless field (jj, the scal- 
ing exponent A is related to the conformal dimension h 
through the relations 



h^h=-A, 
2 ' 

and to the correlation exponent rj by 

1 
A = -n. 
2 ' 



Specific values for 77 relevant to what follows are [g[ rj 

5 

24 



for the Ising model and V — ^ for percolation. 



(2) 



(3) 



In this paper we consider a system at its critical point 
in a domain, with fixed boundary conditions on the 
boundary. Examples of the model system we have in 
mind would be the magnetization m(x, y) of the Ising 
model at its critical temperature defined on the half-plane 
(x,y>0), with the spins on the boundary y — Q kept at the 
value of -|-1, or a percolating system at threshold, where 
m{x, y) is the probability that the site (x, y) is connected 
to the boundary y — 0. We will describe the models in 
this way, even though we will mostly be considering the 
continuum versions of the theories, for which the bound- 
ary condition is that the field is fixed at </> — » cxd at the 
boundary. For the half-plane geometry, the form of the 
magnetization is fixed by translational, rotational, and 
scale invariance 



< m{x,y) >= Cy' 



(4) 



where C is a constant. This result also describes the 
behavior of a lattice model, except when y is as small 
as a few lattice spacings; the different behavior at short 
distances can be viewed as the consequence of corrections 
to scaling, which in any case are missing from Eq. (4). 
Eq. (1) is relevant to a more general problem, because it 
relates the profiles for different-shaped regions; this has 
been used pj to determine the magnetization profile for 
several cases, such as the disk and the infinite strip. 

As will be shown below, conformal invariance implies 
that the critical magnetization profiles for different mod- 
els can be expressed in terms of a single function ?Tii(r): 



TO(r) = mi{r) 



A 



(5) 



where 7Tii(r) is a function that depends on the shape of 
the domain but not on any feature of the model, and 
describes the critical behavior of a mythical model with 
A = I. The meaning of this result is that all confor- 
mally invariant models with the same exponent A have 
the same magnetization profile; then it is enough to look 
at the magnetization profile for the Gaussian models. 
Therefore, Section 2 will begin with a discussion of the 
Gaussian model and end with the proof of the assertion 
(5), and with explicit exhibition of the function mi for 
some simple geometries. 

Section 3 will present some numerical simulations for 
the Ising and percolation models, verifying that they are 
related this way and showing that fitting the numerical 
data to the known form of mi gives an accurate way 
to determine A, in which far from worrying about the 
implications of boundaries and finite size, we are making 
use of their known effects. 



In Section 4 we will demonstrate the usefulness of this 
technique in a determination of the temperature depen- 
dence of rj for the planar spin model. 



II. ORDER PARAMETER PROFILE IN FINITE 

GEOMETRIES ACCORDING TO THE 

GAUSSIAN MODEL 



The Gaussian model has the action 



1 



S^^9 I iV^pdV. 



The order parameter is 



(6) 



(7) 



where the expectation value can be evaluated from the 
functional integral 



m(r) 






(8) 



To calculate this integral we expand the field 9 in terms 
of a set of functions Lpx 



") ^^a\^\{i 



(9) 



and treat the coefficients ax as independent degrees of 
freedom. The boundary condition on spin ordering be- 
comes the condition that ^p\{r) takes a constant value 
(zero, for the moment) for r on the boundary. It proves 
useful to choose these functions to be the normalized 
eigenfunctions Lp\ of the Laplacian 



V^^x{r) = \ifAr), 



(10) 



because these diagonalize the action; the resulting Gaus- 
sian integrals give 



m{r) — e 






(11) 



Interpretation of this result is facilitated by the observa- 
tion that the two-point function 



G(i?,rO==^i^A(i?)¥.A(r) 



(12) 



is the Dirichlet Green function: it satisfies 



V^G(i?,r) = <5(i?,r), 



and gives G{R, f) — for R on the boundary. Thus the 
quantity appearing in the exponent of (11) appears to 
involve limj;^pG{R,r), which does not exist. However, 
this same difficulty occurs in electrostatics in the evalu- 
ation of the electrostatic energy of a unit point charge 



at position f near a conducting boundary: G{R, r) is the 
potential at point R caused by the charge, but to cal- 
culate the energy we must remove the self-energy of the 
point charge. In the present context this suggests that 
we should replace G by 



G{r) = limj,JG{R, r) + —ln\R - r\] 



27r 



and then write 



— c 2q \ > 



m{r) = e ^^ 



(14) 



(15) 



This subtraction can be viewed as a modification of the 
action (6), in which the fields 9 become "normal-ordered" 
[|l[. What has been subtracted in Eq.(14) is independent 
of r - it is an (infinite) constant. The result is that G{'r) 
is finite throughout the domain but diverges to — oo at 
the boundary, thus returning the boundary behavior (5) 
for m{r). We see that by solving an electrostatics prob- 
lem we can obtain the order parameter for the Gaussian 
model defined on the same geometry. 

We now argue that this result describes the critical 
magnetization profile for any model that can be studied 
by conformal field methods. The argument is as follows: 
we can match the order parameter profile for the Gaus- 
sian model in the case of the half-plane geometry (Eq. 5) 
to that of any given critical model by choosing 



1 

Aug 



(16) 



since for this case G{r) — ■^ln\2y\ - the part of the Green 
function due to the "image charge." But then the con- 
formal transformation rule (1) does not depend on any 
property of the ordering field other than the conformal 
dimension h, or equivalently the exponent A. Thus we 
can conclude that any critical model in two dimensions 
will have the same order parameter profile as its match- 
ing Gaussian model, at least for all boundary geometries 
that are related to the half-plane by a conformal trans- 
formation. In view of the simple way that g enters into 
Eq. (15), this result can be restated to say that the or- 
der parameter profile for these cases can be written in 
the form (5). 

Here are some examples: 

*The disk with spins fixed at the value of +1 on the 
boundary: 



<mir)>^C\A--\-^. 



(17) 



Here A is the radius of the disk and r is the distance 
from the center of the disk. This expression is valid |q| 
for both r < A and r > A. 

*The infinite strip r = {x, y) with < y < W , with 
spins fixed at +1 on the boundaries y = and y = W: 



<m{r) >^C{simr^)-^. 



(18) 



The derivation of this expression is discussed in Appendix 
A. 

*The rectangle —^<x<^,0<y< W, with spins 
on aU boundaries fixed at the value +1: 



< m{w) >= C[- 



Im[sn{^w)] 



l-A 



(19) 



where w = (x -\- iy) is a complex variable defined on the 
rectaiigle, and sn, en, dn are the Jacobian elliptic func- 
tions H] . The aspect ratio of the rectangle enters through 
the modulus of the elliptic functions: the quarter periods 
K and K' are related by K' / K = 2W/L. The derivation 
of this result is given in appendix A. Even though elliptic 
functions are a little less well known than the functions 
that appear in Eqs. (17) and (18), the rectangular ge- 
ometry is might be easier to study numerically, and the 
expression (19) is easily evaluated. 

Common to all the problems mentioned above is that 
the boundaries are conformal images of the real axis. 
They can also be obtained M by conformal mapping from 
the half-plane, transforming Eq. (4) by means of Eq. (1). 

The order parameter profile for the Gaussian model on 
an infinite strip with opposed boundaries (spin -1-1 on one 
side, spin -1 on the other) can be found by an extension 
of the method. The Gaussian field 9 is divided into two 
parts: 



e{r)^eo{r) + e,{r), 



(20) 



where Oi{r) is a fluctuating field that vanishes on 
the boundary while 0o(^) = T^y/W is the solution to 
Laplace's equation with the opposed boundaries. Par- 
tial integration then shows 



g /(|V0o|' + |V0i|2)dV, 



(21) 



so that 6i is not coupled to (?o in the statistical averaging. 
It follows that 

m+^{r) ^< i?ee*(^"(*')+^i(^» >== m++(r) cos 7ry/H/, (22) 

where m++ — nii is the profile for the case that the 
boundary conditions are the same on opposite bound- 
aries, which we have already considered. This result has 
been previously obtained [M by means of a conformal 
field theory argument. The argument above now implies 
that a universal prescription can be given for the magne- 
tization profile in a region whose boundary is an image 
of the real axis, with boundary conditions that change 
from -f to -. 

The Gaussian model is known to give the incorrect two- 
point correlation function for the Ising model H , and so 
despite the successes reported above, we should not as- 
sume that the Gaussian model correctly reproduces the 
magnetization profile for all models in all boundary con- 
ditions. A counterexample is the case of the Ising model 



with spins fixed at -1-1 on two ends of a cylinder (i.e. 
a rectangle with periodic boundary conditions joining a 
pair of sides) . Since this geometry cannot be related to 
the half-plane by a conformal transformation (the bound- 
ary is not in any sense connected), the chain of logic that 
allowed us to identify order parameter profiles for the 
Gaussian model with those for nontrivial models is bro- 
ken. The Gaussian model solution for this case is given 
in Appendix A, with the result 

< m{x,y) >^ \e^[^-^y)exp{-^y^T'^, (23) 

where 9i is the theta function |^ of nome q — 
exp{—2TrW/L) and t = 2W/L. As will be seen below, 
simulations do not support this form well, and indicate 
that the Ising model magnetization is better represented 

by 



< m{x, y) >- 



iify)\' 



5(fy)h 



l^i(¥y)n 



(24) 



In comparing these, please note that A = i for the Ising 
model. The logic leading to Eq. (24) is explained in 
Appendix B. 



III. NUMERICAL SIMULATIONS 

We tested the predictions of conformal field theory 
by simulating the Ising model and the bond percolation 
problem on the square lattice. 

For the Ising model we used the Swendsen-Wang 
Monte Carlo algorithm 0] . For each simulation the num- 
ber of Monte Carlo steps per site was 10^. In each simu- 
lation the intial configuration was all spins at the state of 
-|-1. The undetermined constant C was used as a fitting 
parameter. 

The percolation threshold for the bond percolation 
problem on the square lattice is pc — \- We evaluated 
m[r) by averaging 10^ random bond assignments. 

*Disk geometry. We cut a disk from a piece of square 
lattice with unit lattice constant, defined by the sites for 
which \r\ < A. Sites outside this disk were in frozen 
state: all "spin +" for the Ising model or "connected to 
the boundary" for percolation. We averaged the magne- 
tization over all directions to obtain m(r). The results 
are shown in Figure 1, which is a graph of In m(r) versus 
ln{A—r'^ /A). The data is expected to fall on straight lines 
of slope —A, according to Eq. (17). Dashed lines having 
the appropriate slopes have been drawn in for compari- 
son. Since the simulations were done for lattice systems 
and the theory is derived for a continuum model, there 
can be a discrepancy close to the boundary. In our simu- 
lations, the boundary is rough at the lattice scale, and in 
most directions is slightly more than A lattice spacings 
from the center. Therefore we allowed the value used for 



A in the fitting to be slightly larger than the (integer) 
value used to define the disk, chosen so as to eliminate 
curvature that would otherwise appear at the left edge 
of Figure 1 (order close to the disk boundary) . We found 
the optimal value to be ^ = radius + 0.7 for the Ising 
model and A = radius + 0.5 for percolation. By fitting 
the data for the disks of three different radii and finding 
the average slope we get A — 0.1234 ± 0.0005 for the 
Ising model, and A — 0.1036 ± 0.0002 for percolation. 
The theoretical prediction from Eq. (17) is A = 0.125 
for the Ising model, and A « 0.10417 for percolation. 
Thus this seems to be a viable way to determine A. 
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FIG. 2. Average value of magnetization as a function of 
the distance from the boundary for L=W=100. 



The dashed lines arc the theoretical curve 



< m{r) >~ C{simT — ) «, 



(25) 



obtained from Eq. (18), displaced by a small con- 
stant amount upwards and downwards, since otherwise 
the theory and data would greatly overlap. The lower 
part of this figure is a separate graph of the same in- 
formation, now presented as the relative error ER = 
100 * {rrith — 'rnexp)/'mexp- We will use this format in 
most of our figures. 

*Ising model on infinite strip with H — boundaries. As 
above, we approximate the infinite strips with rectangles 
(0< x< L,0< y < W),L>W. The results for the case 
L=W=100 arc shown in Fig. 3. 



ln(A-r VA) 

FIG. 1. Linear fits for the Ising model and percolation on 
the disk 



*Ising model on infinite periodic strip with -I— I- bound- 
aries. In a simulation we can work only with finite strips. 
A rectangle (0< x< L,0< y < W) which is periodic in 
X direction approximates an infinite strip really well if 
the periodic dimension (in our case L) is greater or equal 
then the dimension in the nonpcriodic direction (in our 
example W). The results for the case L = W are shown 
in Fig. 2. The spins on the boundary y=0 and y=:W 
were fixed at the value of +1. 
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FIG. 3. Average value of magnetization as a function of 
the distance from the positive boundary for L=W=100. 

The spins on the boundary y=0 were kept at the value 
of +1, while the spins on the boundary y=W were fixed 
at the value of -1. The dashed lines are the theoretical 
curve: 



V —1 y 

< m{r) >— C(sin7r — ) ** costt 



W 



w 



(26) 



obtained from Eq. (22), displaced as before. 

^Percolation on infinite periodic strip. As in the Ising 
model, we approximate the infinite strip with a square 
(0< x< 100, 0< y < 100) which is periodic in x direction. 
The results are shown in Fig. 4. 
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FIG. 4. Order parameter as a function of the distance from 
the boundary for a periodic square 100x100. 



The dashed lines are the theoretical curve: 

< mir) >= C(sin7r-V7)"^- 
14/ 



(27) 



*Ising model on rectangle. The rectangle (0< x< L,0< 
y < W) was defined by spins fixed at the value of +1 on 
all boundaries. The average magnetization was measured 
as a function of the distance from the boundary y = 
O.The results for the square (L,W)=(100,100) are shown 
in Fig. 5 and in Fig. 6. The magnetization in Fig. 5 was 
measured along the line (x = 10, y) while in Fig. 6 the 
measurement was along the line (x ~ 50, y). 
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FIG. 5. Average value of magnetization along the line 

(x=10,y) 
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FIG. 6. Average value of magnetization along the line 

(x=50,y) 



The dashed lines are the theoretical curve: 
< m{w) >= C[- 



Im[sn{^w)] 



-\cn{^w)dn{i^w)l 



(28) 



obtained from Eq. (19). 

*Percolation on rectangle. The rectangle is (0< x< 
100, 0< y < 100). The order parameter is given as a 



function of the distance from the boudary y = 0. The 
results along the line (x = 10,y) are given in Fig. 7 while 
in Fig. 8 the measurement was along the line {x ~ 50, y). 
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FIG. 7. Order parameter along the line (x=10,y) 
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FIG. 8. Order parameter along the line (x=50,y) 
The dashed lines arc the theoretical curve: 
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(29) 



*Ising model on the periodic rectangle. We used the 
rectangle (0< x< L,0< y < W) periodic in the x di- 
rection, with the spins kept at the value of +1 on the 
boundaries y=0 and y=W. As we have seen above, the 
case L>W resembles the infinite periodic strip. Here we 
will consider L<W. In Fig. 9 we show the average value 
of magnetization for rectangle L=80,W=101. In Fig. 10 
the rectangle is L==40,W==101. 
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FIG. 9. Average value of magnetization for the (80,101) 
rectangle. The data is displayed twice, with a vertical dis- 
placement. The upper set of lines compares the magnetiza- 
tion to the prediction (23) based on the Gaussian model. The 
lower set is Eq. (24), supplied by conformal field theory. 
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FIG. 10. Average value of magnetization for the (40,101) 
rectangle. Same comparisons as in figure 9. 

In the lower part of the figure, the full line is the rela- 
tive error for the conformal field theory prediction while 
the Gaussian model error is represented by the dashed 
line. The same convention will be used in Fig. 11 and in 
Fig. 12. The agreement between the data and the pre- 
dictions of the Gaussian model is not very good; Eq. (24) 



coining from conformal field theory does much better. 

*Percolation on the periodic rectangle. The position of 
percolation theory in conformal field theory is as yet un- 
clear: it apparently is not a unitary minimal model, and 
it is not known whether there is a differential equation 
to which TO(r) is a solution. In view of the failure of the 
Gaussian model to describe the magnetization profile of 
the Ising model, we have no theory for this case. We have 
chosen to give the same analysis as was done for the Ising 
model, changing only the exponent A. Thus the top set 
of curves in Figures 11 and 12 compare the behavior of 
the percolation order parameter to the Gaussian model 
(23), while the lower set compare the same data to 



< m{x,yj >= ( -— j-^ )6, 



(30) 



which is related to the Ising model expression (24) 
through the universality hypothesis Eq. (5). There's no 
reason for this to work! But it seems to agree with the 
data! 
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FIG. 11. Order parameter for the (81,101) rectangle 
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FIG. 12. Order parameter for the (41,101) rectangle 

IV. PLANAR SPIN MODEL 

The planar spin model consists of unit spins that can 
point in any direction in the x-y plane. The Hamiltonian 
for this model can be written as 



H = -J Y^ cos{6, - Oj), 



(31) 



<i,3> 



where the sum is over the nearest neighbors, di is the 
angle the i th spin makes with the x axis. 

The model undergoes a Kosterlitz-Thouless (KT) tran- 
sition |jl3-n3] at the critical temperature T — T^. Below 
this temperature there is algebraic order at all temper- 
atures, implying a lack of a length scale and thus con- 
formal invariance; unlike three dimensional phase transi- 
tions, however, the exponent rj depends on temperature. 
We determined m(f) by means of simulations on a disk 
of radius 100, with 9 ~ Q on the boundary (i.e. the 
spins on the boundary are constrained to point in the x 
direction), using Wolff's algorithm jl^l. After the equi- 
libration we performed 10® Monte Carlo steps per site. 
By using < m(r) > obtained from numerical simulations 
we can find rj, since Eqs.(3) and (17) imply 



In < m(r) >- 



^iMA-^) 



constant. 



(32) 



Linear regression was used to determine the slopes of 
the lines in Fig. 13, which gives the exponent ri{T). 




FIG. 13. Linear fits for temperatures 0.85J, 0.80J, 0.75J, 
0.70J, 0.65J, 0.60J, 0.55J, 0.50J 

Our results are given in Table |: 



T 0.85J 0.80J 0.75J 0.70J 0.65J 0.60J 0.55J 0.50J 



rj 0.2187 0.1882 0.1660 0.1521 0.1348 0.1201 0.1103 0.0951 



TABLE I. Critical exponent as a function of temperature 

According to Kosterlitz-Tliouless theory ll^, the long 
wavelength spin fluctuations are described by the Haniil- 
tonian (32) with a renornialized coupling constant J(T), 
given by 



J{T)/T ■ 



(constant) 



T-Tr. 



Tc 



(33) 



where the constant is determined by the core energy for 
a vortex. This determines the temperature dependence 
of 7] — T/2ttJ, so that close to the Kosterlitz-Thouless 
transition 



[\-v{T)r 



Tc-T 
Tc ' 



(34) 



implying that rj approaches its critical value by a square- 
root cusp. We note that many of the previous studies 
[ p5| , p8[ of 77 (T) have incorrectly assumed a linear depen- 
dence. 
Figure 14 shows the analysis of the data from Table |. 
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FIG. 14. Data comparison 

The points follow a straight line as predicted by 
eq.(34). The critical temperature is given by the intersec- 
tion of this line with the temperature axis and is found 
to be 



rc/J = 0.863 ±0.009. 



(35) 



In Fig. 14 we also plotted the results obtained by To- 
bochnik and Chester (TC) [|l5|. We notice that their 
data consistently fall below our line, but are consistent 
with it. TC analyzed their data using a linear extrap- 
olation to obtain Tc/ J = 0.89 ± 0.01. Our reanalysis 
of their data indicates a lower value. We also quote the 
estimates for Tc obtained by some other authors: We- 
ber and Minnhagen |l|l give Tc / J = 0.887 ±0.002 while 
Olsson II finds Tc/J = 0.8922(2). 

The results obtained by computer simulations in the 
critical region are affected by the finite size effects - one 
is trying to obtain predictions for the infinite system by 
simulating a finite system. Our approach avoids this 
problem because the theory says that the critical expo- 
nent rj for the model defined on a disk must in principle 
be equal to the exponent 77 for the infinite system. The 
only approximation involved is the passage from a dis- 
crete system to continuum. 



V. SUMMARY 

In this work, we studied the magnetization of the Ising 
model and percolation defined on finite geometries. We 
derived the order parameter for rectangle by using con- 
formal symmetry. We showed how the results previously 
obtained by using conformal field theory can be obtained 
by studying a simple Gaussian model. We numerically 
studied order parameter in various geometries by using 
Monte Carlo methods and compared analytical results 
with the outcome of our simulations. 



VI. APPENDIX A 

To determine the order parameter for the model de- 
fined on an infinite strip ( x, 0<y<W) with spins fixed 
at the value of -1-1 on the boundaries, we find the electro- 
statics Green function for the infinite strip with Dirichlet 
boundary conditions: 



G(Z,.) = ^(ln|sm^ 



— In I sin 



ttIZ- 



2W 



-1), (36) 



where z = x + iy is the complex representation of the 
plane, and z* — x — iy is the position of an image of z 
p2[ . This result can be obtained by conformal mapping 
from the half-plane; or from the observation that Eq.(36) 
has logarithmic singularities of opposite signs at the po- 
sitions of all images (in the two boundaries) of z and is 
a solution to the two-dimensional Laplace problem else- 
where, since it can be written as the sum of a function 
of Z and a function of Z* . The non-singular part of the 
Green function is then 



G{r 



1 , , . Try 



(37) 



Eq. (15) then gives the result quoted in text. 

The magnetization profile for a rectangle can be found 
using a conformal mapping, as follows: The semi-infinite 
complex plane z = (x, y > 0) is mapped onto the rect- 
angle (—-J < u < ^,0 < iv < iW) in the complex plane 
w = (w, iv) by the Schwarz-Christoffel transformation 



wiz) = -— / g{z')dz', 



where 



with 
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\g{x')\dx' , 



(38) 



(39) 



(40) 



K' = 



\g{x')\dx'. 



(41) 



The ratio of K and K' is related to the nome q of the 
Jacobian elliptic functions |^ 

q^e-''^. (42) 

The inverse of the function w{z) is the Jacobian elliptic 
function z = sn{^w), which implicitly depends ||6| on 
the nome q. Then we can put (4) and (38) into Eq. (1), 
'"^^-w)], to obtain 



noting that y = Im[sn{- j^ 
< m{w) >= C\- 



Im[sn{^^w) 



1 _ k^sn^^w)Jl - sn'^i^w 



We use the relations between Jacobi elliptic functions sn, 
en, dn H] to obtain Eq. (19). 

This result can also be obtained through the connec- 
tion with the Gaussian model. This requires finding the 
Dirichlet Green function for the rectangle, which can be 
found by the method of images. The images of the point 
r lie on four rectangular Bravais lattices, and the Green 
function itself is an infinite sum of the corresponding log- 
arithmic terms. The sum can be evaluated by observing 
that exp(2TrG{R, f)) has only simple poles and zeroes, 
which form periodic arrays, and thus it is a function of 
elliptic type M. Inspection then gives 



exp{2TrGiR,r)) 



e,{n^)e,{n^ 



■) 



^i^^Wii^^^y 



(44) 



2L )"i-\" 2L 

where Oi is the thcta function |1^,^, Z is the complex 
representation of R, and the z's are the complex repre- 
sentations of r and the image points z' = z* , z" = L — z, 
z'" — L — z* . The resulting expression for the order 
parameter profile is 



< m{x,y) >= I 



02(^^) ' 



(45) 



This representation is equivalent to the one used above, 
but is easier to evaluate numerically, since the series 
representations H for the theta functions converge ex- 
tremely rapidly. 

The Dirichlet Green function for a i x M^ rectangle 
with fixed boundary conditions on two sides and peri- 
odically connected along the other two has logarithmic 
singularities at Z = z + mL -\- inW, and logarithmic sin- 
gularities of opposite sign at Z = z* + mL + inW, where 
z* — X — iy. This behavior is reproduced by the function 



G{Z, z) 



1 



-ln\ 
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(46) 



Here t = 2W/L specifies the shape of the rectangle. We 
take the limit Z i-^ z as we did in the case of nonperiodic 
rectangle and obtain eq. (23). 



(43) 



VII. APPENDIX B 

According to conformal field theory pi, the N-point 
correlation functions of many critical models can be 
found by solving a differential equation originating in the 
algebra of the operators corresponding to the fields. In 
many cases, the correlation function factors into a prod- 
uct of holomorphic (depending on z) and antiholomor- 
phic (z*) parts, or a sum of a small number of such 
terms. Cardy Q has given a prescription for introducing 
a straight homogeneous boundary into this formalism. In 
effect, the N-point function is found by solving the 2N- 
point differential equation for the holomorphic part, and 
then evaluating this with the second set of points at the 
image position (i.e. zjv+i = z^ if the boundary is the real 
axis). Since the differential equations are of higher than 
first order, there are several solutions to choose from, 
and the choice is determined by the boundary condition 
(e.g.Dirichlet) desired. 

The problem we wish to consider is the rectangle with 
Dirichlet boundary conditions on two opposite edges and 
periodic continuation on the other two edges. This 
is equivalent to a finite length cylinder with Dirichlet 
boundary conditions on the ends, or to a torus that is 
cut by a single Dirichlet seam. 

The differential equation appropriate for the Ising 
model spin operator on a torus has been given by Eguchi 
and Ooguri pi| , and the solutions for the two-point func- 
tion are given by Di Francesco et al. |22]. They all have 
the form 



< aa >= 



[Oi{zi - Z2)]5 



(47) 



where 9i, [v = 2,3,4) are the theta functions with the 
periodicity of the torus. The combination < aa >2 + < 
aa >3 has the desired periodicity in Im z. We then follow 
Cardy by evaluating this with zi = z, Z2 — z* . 
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